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Abstract 

From a refined non-relativistic quark model that fits the baryonic low- 
energy spectrum the study of strong pion decay processes within an ele- 
mentary emission model scheme points out the need of incorporating size- 
contributing components into the baryon wave functions. In particular the 
effect of a (qqq qq) component is investigated in the framework of a quark 
pair creation model. 



Preprint FTUV/95-71, IFIC/95-74 
Preprint |arch-ive/ 9606038 



1 Introduction 



In QCD, the basic theory of strong interactions (responsible for the hadronic 
structure), elementary constituents (quarks) which are massless if chiral sym- 
metry is preserved can acquire a dynamically generated mass associated to the 
spontaneous breakdown of chiral symmetry that takes place when the strong 
interaction reaches a critical intensity. On the other hand, although not conclu- 
sively proved in (3+1) dimensions, QCD seems to be a confining theory. Assum- 
ing that the energy scale for the chiral symmetry breaking is higher than the 
corresponding one for confinement (~ Kqcd)-, one gets the image of a baryon as 
a confined system of interacting massive quarks and Goldstone bosons coming 
from the symmetry breakdown. 

Little more than this qualitative statement can be said from QCD at the 
current moment since the lack of an adequate method to describe a bound state 
within a quantum field theory, together with the non-perturbative character of 
the interaction below Aqcd, makes difficult any attempt to go further. As a 
consequence several models of hadronic structure have been developed. Among 
them perhaps the ones closer related to the image of QCD just drawn are the 
bag models where the quarks move almost freely inside a confinement region, the 
bag, surrounded by a mesonic cloud. However the separation of the center of 
mass motion (indispensable to tackle the problem of baryon structure) added to 
the fact that one has to simultaneously solve for the movement of the quarks and 
the bag surface, make the model technically complex. To avoid these difficulties 
one could wonder whether it would be possible or not to construct an equiva- 
lent (regarding the low-energy description) non-relativist ic model. Undoubtedly 
the extensions and refinements from the naive quark model of hadron structure 
[|TJ are steps in this direction. In this line De Riijula et al. [0 derived an in- 
terquark potential from the non-relativistic reduction of the one gluon exchange 
diagram in QCD to be considered together with the phenomenological confining 
interaction. The resulting non-relativistic quark model (NRQM) reproduces the 
image of the quark core of the bag, the mesonic cloud having been obviated. 
The absence of mesons together with the non-relativistic treatment poses serious 
doubts about its usefulness unless the incorporation of new terms in the poten- 
tial and/or the effectiveness of the parameters of the model allows to mitigate at 
least partially these shortcomings. As an immediate consequence of this way of 
proceeding, the constituent quarks loose a direct connection to the quark fields 
of the QCD lagrangian, representing effective degrees of freedom. This viewpoint 
is also supported by the absence of a strong necessity to introduce in the interac- 
tion between quarks some terms issued from the one gluon exchange such as the 
spin-orbit or tensor forces |3|, . 

The first requirement to such an effective model should be a good description 
of the hadronic spectrum (at least at low energies) that allows a correct assign- 
ment of the theoretical states to the experimental ones. This can be accomplished 
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with relative success for mesons as well as for baryons by a "minimal" model in- 
volving confinement plus part of the one gluon exchange (OGE) potentials ||]. 
Nevertheless an endemic problem in the low energy spectrum remains, that is 
the impossibility to correctly predict the masses for the first radial excitations 
in the baryonic octet and decuplet. The systematic of this wrong prediction led 
some authors to add a phenomenological three-quark force getting the cor- 
rect baryonic spectrum (for negative and positive parity states with and without 
strangeness simultaneously) up to an excitation energy of 0.7 GeV without need 
to resort to any specific resonance mechanism. This model, that will serve as our 
starting point now, presents however some conceptual problems which we shall 
consider firstly. 

The second requirement is the correct description of the baryon properties, 
what involves the wave functions of the baryons on the one hand, the operator 
relative to the process under discussion on the other hand. Having a model 
that correctly describes the low energy baryonic spectrum, it is appropriate to 
compare its predictions to those of other models that do not do so well for the 
baryon masses. This may allow us to determine whether there is some bias to be 
introduced in the comparison to experiments due to a poor prediction of baryon 
masses (though some discrepancy may be due to different physical ingredients). 
This comparison was our main intent originally for looking at the strong pion 
decays of baryons. It however appeared that results were very much sensitive to 
the microscopic description of the process. This led us to examine various aspects 
of the problem, where relativity is often present in one way or another. The 
present work is part of a general program tending to extract from the comparison 
of the models to each other and to the data some general features of the "true" 
wave functions that in their turn should serve as a guide to further refinements of 
the models themselves. The contents of the presentation are organized hereforth 
as follows. In section 2 we review the quark models we make use of and their 
predictions for static properties. In sections 3 and 4 we study strong pion decay 
processes within two different approaches, namely the elementary emission model 
(EEM) and the quark pair creation model (QPCM) that we modify in order 
to correct its energy dependence. Finally in section 5 we summarize our main 
conclusions 

2 The Quark Models 

The non-relativistic quark model of hadron structure describes the hadron (color 
singlet bound state) in terms of a definite number of components (constituent 
quarks q, and antiquarks q) that interact through an effective potential. The 
analysis of the mesonic (qq) and baryonic (qqq) spectrum leads to a minimum 
two-body potential containing the basic QCD motivated, confining, coulombic 
(OGE) and spin-spin (OGE) qq interactions, of the form: 
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where D is a constant to fix the origin of the potential, is the distance between 
quarks i and j and a denote the Pauli matrices. The Yukawa form of the spin- 
spin term replace the 5(f) contact interaction of the OGE potential to avoid 
an unbounded spectrum when solving the Schrodinger equation [Q. The mass 
parameters m^j are chosen to fit the baryon magnetic moments and the three 
parameters a 2 , k = K a , r , fitted from the meson spectrum J|] (see table 1) 
provide a good description of the baryonic octet and decuplet || , a part of which 
is reproduced in fig. 1. A look at this figure makes clear a general deficiency 
in the description: the masses of the positive parity excitations that correspond 
to radial excitations of the totally symmetric spatial components of the wave 
function are systematically higher (from 200 to 400 MeV) than the experimental 
ones. Any attempt to correct this situation by a refitting of the parameters causes 
undesirable effects on the rest of the spectrum. Then the open question refers 
to the possibility of solving the problem by means of the introduction of some 
physically founded new term in the potential. 

Attending the mixed symmetric spatial structure of the first negative parity 
excitations it becomes obvious that they will be much less affected by an inter- 
action tending to group the three quarks than the spatially symmetric ground 
states and its first radial excitations. Having in mind a genuine three quark in- 
teraction as the exchange of two sigma mesons at the same point, Desplanques 
et al. J5| proposed a three-body potential 

1 TA p-monj p -m r ik 

Vf = - V Vo (2) 

2 i ^ k ^ i m i m j m k m Q r i:j m Q r ik 

The low energy baryonic spectrum can be then well reproduced with a 'two 
plus three body' interaction: 

Vn = V}? + V;f (3) 

(2) 

where Vjj differs from Vj in the values of the parameters. As a consequence the 
unified description of both the mesons and baryonic spectra given by Vj is lost. 
More than a drawback one should consider this as the natural outcome of the 
effectiveness of the parameters (through them for instance, some mesonic effects, 
only present in the baryonic case might be incorporated). 

The form used in ref. [[J for the potential, eq. (|3|), shows some non-appealing 

(2) 

features. Concerning the two-body piece v;; ; , the range of the spin-spin interac- 
tion, r = 1.28 fm, needed in || is very far from the one which could reasonably be 
expected given the original 8(f) form of this interaction 0. This comes from the 
minimum number of parameters' constraint k = K a . Nevertheless one can get a 
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much more "reasonable" range (0.4 fm) without increasing the number of param- 
eters by recovering the original OGE potential relation between the coulombian 
intensity k (k — — |a s ) and the spin-spin one n a = — |a s ), i.e. K a = |/c. The 
resulting value of tq might represent somehow an effective average between fm, 
the OGE 5-term range, and ~ 1 fm, the range of the quark pionic exchange not 
implemented in our model. This impression seems to be reinforced by the corre- 
sponding value of the effective QCD coupling constant, a s ~ 1, that doubles the 
one needed in models where pionic quark exchanges are explicitly incorporated 

i- 

From now on we shall call Vjj the two plus three body interaction with Vj\ 
corrected as explained. The fit to the spectrum with this corrected potential (see 
table 1) a part of which is reproduced in fig. 2, hardly differs from that in ref ||. 

Regarding the 3-body piece of the interaction (H), Vjj, its long range, 0.8 fm, 
can find a justification attending the small value of the coefficient of the linear 
potential (or equivalently its small contribution to the string tension). For higher 
angular momentum states where the quarks are very far apart, the long range of 

(3) 

Vjj would allow a contribution to the string tension to explain the leading Regge 
trajectories. It is difficult however to justify if we think of the three nucleon force 
this 3-quark term would give rise to. To correct for this behaviour a less singular 
form of the interaction is required. For the sake of technical simplicity we shall 
adopt here a gaussian form 

V® = V expl-X:^l (4) 

very easy to be used in the hyperspherical harmonic formalism (Appendix A) M 
we work with. 

The parameters of the potential 

Vni = Vifj + Vjfj (5) 

appear in table 1 (note that the range of V® has gone down to 0.25 fm and 
that we have relaxed the constraint between k and K a ), and the comparative 
spectrum, of similar quality to the V/j's one, appears (in part) in fig. 3. 

With regard to other parameters of the potential some comments are in order. 
For the quark masses, m u , m<j, a value between 300 and 360 MeV, that fits 
the nucleon magnetic moments (assuming an SU(6) spatially symmetric wave 
function) within a 10 % of error, is generally accepted. It can be shown (Appendix 
A) that, from the potential we are dealing with, one can obtain via an adequate 
redefinition of the parameters of the potentials, a non-strange energy spectrum 
(for quark mass m') whose energies, E', are related to the initial ones E (for 
quark mass m) by 
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E' = —E (6) 
m' 

and with no other change than a distance rescaling in the wave functions. Hence 
given a fit to the spectrum for a given quark mass, eq. @ gives the shifted relative 
energies when changing the quark mass, i.e. the baryon magnetic moment. 

The strength of the confining term is considerably reduced due to the presence 
of the three body force (see table 1). This gives rise to a grouping of the states in 
the 'high' energy part of the spectrum that could be suggesting the need of some 
qualitative change in the three-body or in the confinement pieces. Nevertheless 
the lack of precise experimental information in this region together with the very 
probable presence of significant relativistic effects prevent us from extracting 
quantitative conclusions from it. 

Once, via the fit of the spectrum, some physical states are unambiguously 
ascribed to model states (up to 0.7 GeV excitation energy) we should focus on the 
description of the baryonic properties. As explained before the nucleon magnetic 
moments are reasonably well reproduced (see table 2) with masses between 300 
and 360 MeV indicating the dominance of the spatially symmetric component of 
the wave function (see fig. 4). A different case is the mass and charge root mean 
square radius (r.m.s.). For the nucleon charge radius for instance the predicted 
values are too low compared to the data. This can be understood by taking into 
account that our effective scheme only incorporates the description of the core of 
quarks, whose pretty small size has to do with the pressure exerted by the mesonic 
cloud that in our meson-absent models is simulated through the effective values 
of the parameters. Then, if the mesonic contribution to the total mass is small 
as compared to the core contribution (as it happens to be the case in the bag 
models), one can understand the r.m.s. -spectrum puzzle, i.e. the impossibility 
within a model scheme that only plays with effective quark degrees of freedom of 
simultaneously fitting the spectrum and the total sizes. 

Furthermore the absence of (qqq qq) components in the wave function should 
be at least in part also responsible for the high values (5/3) predicted for the 
axial coupling constant qa as it is suggested by bag model calculations |L0| . 



As a corollary it may be established that when dealing with physical pro- 
cesses where the mesonic cloud plays a relevant role as such the use of a spectro- 
scopic quark model requires the use of effective operators that can give account 
of it. Hence for the strong pion decay we consider next a careful analysis of 
the transition matrix elements is needed before a comparison between different 
spectroscopic quark models makes any sense. 

3 The Elementary Emission Model (EEM) 

To study strong baryonic decays we shall follow the elementary emission model 



approach developed long time ago Jll|| . The decay takes place through the emis- 
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sion of a point-like pion by one of the quarks of the baryon (fig. 5) and therefore 
the coupling constant is unique for all the processes. 

Up to order (^ SL ) the matrix elements for the process B — > B'tt can be ex- 
pressed as (B'\H\B) fg with 



u = 3i 1 fwar a (3)u 

(2tt) 3 / 2 (2^)V2 ' 



a^ke-^ - ^- {a ^ 3 \ e -^} 
2 Trip 



(7) 



where the factor 3 comes from the number of quarks, u n and k stand for the 
pion energy and trimomentum respectively and r 3 is the quark 3 coordinate. The 
upper index (3) on the operators of spin, a, isospin, r, and quark trimomentum, 
p, indicates that we have chosen for later technical simplicity the quark 3 to 
perform the calculations. The coupling constant f qqn is related to the usual g qqTC 
through the quark-pion mass ratio: 

m n 2m q 

and is the only free parameter of the model. 

The 'non-relativistic' qqir interaction given by H may be obtained from the 
qqn invariant interaction with the pseudovector coupling C qq7T ~ ^> q f^^^! q d^n. 
It differs from that obtained with the pseudoscalar coupling C qq7T ~ ^fgrj^^qTr 
since at first order in (—), this one does not give rise to the term proportional to 
uj n in eq. (0) (it only gives the ak term). The examination of the higher order 
terms in (^) in the pseudoscalar coupling case shows that there are some (|^) 2 
corrections to the ak term which are identical to those that would be derived 
in the pseudovector coupling. As an expansion in (— ) 2 is meaningless in view 
of the large value taken by this quantity (larger than 1, with often a destructive 
interference at the first order) we prefer to freeze these p 2 terms and replace them 
by some constant. In practice, this means that they will be hidden in the coupling 
fqqn to be fitted to reproduce the /nn-k coupling. The other (^-) 2 corrections in 
the pseudoscalar coupling represent the difference in the kinetic energies of quarks 
in the initial and final states, which enters in u n (u^ = Eb — Eb>)- The difference 
in potential energies in the initial and final states, which is also included in u n , 
is provided by terms involving the excitation of a quark-antiquark pair (this one 
is known to be important in the pseudoscalar coupling). 

Notwithstanding that it corresponds to a higher order term in (^-) in the 
pseudoscalar coupling, we shall maintain the term proportional to u n in eq. (0), 
because of its different functional structure. Altogether, eq. (|7|) amounts to retain 
for the different types of terms, ak and ap , corresponding to the lowest order 
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contributions in (—), relativistic corrections being embedded in the parameter 
f qqn or in the factor u n itself. 

In order to evaluate the matrix elements of the hamiltonian H the wave func- 
tion of a baryon of spin J and third component J z is factorized as a plane wave 
of trimomentum P (the baryon trimomentum) for the center of mass motion 
multiplied by the intrinsic wave function, i.e. 

\ B ) = TtA^^ MS, &)®b(S, M s - /, Mj)] i j (9) 



(2tt) 3 / 2 



JJz 



where stands for the spin-isospin part and \1/_b is the spatial internal wave 
function in terms of the Jacobi coordinates £1, £2- 
Then the amplitudes for the decay, defined by: 

J/>X (B'\H\B) JA = jA^KPb - Pb> - k)4 J) (10) 

where % and / refer to the internal quantum numbers of the initial and final 
baryons, are given by: 



A 



3^ /, 



qqir 



( T <* (3)\t 



J'X 



V 6m„ / 



2m„ 



^b^i,^b(S,Ms-J,Mj 



,/a 



(11) 



The dependence on Pb may suggest a dependence of the amplitudes on the 
frame of reference. In the present case, the corresponding term combines with the 
first term proportional to k in eq. ( |TTj) to get a term proportional to k — -^-Pb 

(~ k — m^j^) in the non relativistic limit), which ensures the Galilean invariance 
of the amplitude at the lowest order in This invariance of the amplitude at 
the quark level was imposed by Mitra and Ross [HJ to introduce in the qqw 
interaction some u n dependent term, quite similar to the one in eq. (|7|). We shall 
work in the frame where the decaying resonance is at rest (Pb = 0). For later 
purpose it is convenient to reexpress the amplitudes in the more compact form: 



Ay 



3^ fqqir 

(2^)V2^ 



1 + 



6m„ 



N, 



i,f 



2m n 



-R 



(12) 



The notation Nij, Rij is a technical one. R includes derivatives (coming 
from the momentum operator p^ 2 ) whereas N does not. They are directly related 



S 



to the more physical recoil (depending on the recoil momentum (p^ 3 ^ — k)) and 
direct (depending of the pion momentum k) contributions. 

Finally the decay width is obtained from the amplitudes in the initial baryon 
(B) rest frame as ||12|| : 

where the overline indicates the average over the initial spin-isospin. The kine- 
matical magnitudes and the baryon masses depend on the spectroscopic model, 
so that for Vu and Vm they are pretty close to the experimental ones, whereas 
for Vj the difference (especially for the Roper decays) can be very significant. 
Nevertheless as we shall be interested (once the transition operator has been cho- 
sen) in testing the baryonic wave functions we shall use the same kinematics, i.e. 
the experimental values of tub, Eb>, u n and k, in all our calculations. 



3.1 Fitting of /, 



qqir 



There is no unique criterium in the literature to fix f qqn . For instance f qqn could 
be chosen to get the best global fit to the data JHJ . Nonetheless, for the purpose 



of analyzing the baryonic wave function another possibility is more convenient, 
namely to fit f qqn to reproduce the NNir interaction at low momentum transfers. 
Thus we compare the matrix element for the process p |— > p j 7r° calculated at 
the quark level from H, eq. (0), with the same matrix element at baryonic level 
obtained from: 



H 



Bar 



fNNiv{k) 



(2tt) 3 / 2 (2w 7r )V2 m , 



[cT N k)e 



-ikR 



2m N L > 



(14) 



By assuming the nucleon to be a totally symmetric spin-isospin SU(6) state 
(the mixed-symmetric component of the nucleon has less than 2 % probability 
for our quark models) one gets: 



fNNTr(k) — fc 



5 (1 + ^ 

3 (l + JihL. 



F(k) 



(15) 



where F(k = \k\) contains all the information about the spatial structure of the 
nucleon as provided by the quark model. Explicitly: 



(16) 
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where X represents an integral of the general form: 



= / « ^(0^(^0^.(0 (17) 

J TO standing for the Bessel function of order m and (^c ) f° r the hyperradial 
part of the wave function for the i (j) channel of baryon B (C). 

In the limit (k — > 0, u n = m n ), the spatial structure reduces to F{k) — > 1, 
and then 

, _ 3 fNNTrjO) n „x 

J ^ ~ 5 B(m w ) [ ' 

having defined 

= ^— ^ (19) 

Hence the only dependence of f qq7T on the spectroscopic model comes from the 
non-relevant quark mass differences, having the values f qqn = 0.602, 0.604, 0.600 
for Vi, Vn and Vm respectively (/awv = 0.998). 

Once f qq7T has been fitted in this almost model independent manner, the pre- 
dictions of the models serve as a test, within a EEM scheme, of the baryonic wave 
functions. For instance, by writing Jnntt in terms of a normalized form factor 
G(Jfe): 

fNNn(k) = W(0)G(fc) (20) 
with G(k = 0, u n = m n ) = 1, we have 

G < fc > = F «|Sy < 21 > 

The results appear in fig. 6, against the phenomenological parameterization 
of ref. ||15|| . As can be seen neither Vu nor Vni provide with a good form factor 
(although the deviation is at most of a 10 % for small k values). The reason for 
it is clear when making an expansion in powers of k: 

G{k) = l-\k\r 2 ) + ... (22) 
6 

Thus for the small values of k one is directly testing the root mean square 
mass radius: 



, 2) = ( Eii^(^__R)2 = i 
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which turns out to be too small for Vji and Vj-jj. This is not certainly a surprise 
from our considerations in section 2 since it can be an indication that relativistic 
effects, pionic components (qqq qq), etc may be playing a role. Alternatively it 
may be suggesting the need to refine the three-body potential. 

3.2 Strong pion decays 

To pursue the analysis of the wave functions, within the EEM, we evaluate the 
widths for strong pion decays of nucleon and delta resonances. The results are 
listed in table 3. Although it is difficult to extract general features from them (the 
almost zero widths for the processes involving the ground state and its first radial 
excitation, A/" (1440) — > Ntt and A(1600) — > Air, comes from the orthogonality 
of their radial wave functions) we can say that for decays clearly dominated by 
the non-derivative part of the amplitude (JV(1520) -> Ntt, A(1232) -> Ntt) the 
three potentials work reasonably well, the better the result the bigger the pre- 
dicted radius for the nucleon. This comes along the same line that our previous 
discussion for the NNtt form factor since the non-derivative terms involve the 
average wave function overlap. In order to try to include size contributing com- 
ponents, we shall explore next the possibility of incorporating the pion structure 
in a non-relativistic scheme. 



4 The Quark Pair Creation Model (QPCM) 



To implement the mesonic structure several formalisms have been developed [fL2 
all of them sharing the image of the meson emission as the creation of a qq pair 
that by later recombination gives rise to the outgoing meson. 

At the effective level we work, the choice of one or other QPCM comes mainly 
motivated by simplicity. In this sense the 3 Po QPCM is manageable and compa- 
rable in the limit of point pion to the EEM. Schematically the process is pictured 
in fig. 7. The qq pair created has the quantum numbers of the vacuum: fla- 
vor and color singlet, zero momentum and total angular momentum J PC = ++ 
(=>• L — 1, S — 1). This translates into a transition operator: 



T = - S / d Pq d Pq 

i,3 



3l5(p q + p q ) £(110|m, -m)yr(p q - p q )Zr. 



bl(p q )d)(p q ) (24) 

where 7 is the (dimensionless) parameter strength of the model, yff is a solid 
harmonic and Z~™ contains the color-spin-isospin wave function of the pair. 
The matrix element for the process B — > B'M is then written as: 
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(B'M\T\B) = -3 1 Y.( 110 \ m >- m )(®Bi$M\$B<S>;™)I m (B-,B'M) (25) 

m 

where $ stands for the spin-isospin wave function and 



dp x dp 2 dp 3 dp 4 dp 5 y^(^ A - ^ v )5(^a + " v ) ' 

V*b> (Pi. P2, P^liiPs, Ps)^b(pi,P2, Pa) 
5{P B - P B , - k) [ d£ dg 2 d£&%,(g 1 ,&0& 1 &)* B (£ 1 ,£&6) 



(2tt) 3 / 2 

the kernel 0(^2, £2) being the non-local operator: 



aw-y ffo +?%) - ^+ j^B) *ir(y f (S - &) 

exp \i\l\p B {& -6)) exp Ulfc(g -g) j exp ) (27) 



in terms of the meson wave function ^ m- 

Let us note that, in the rest frame of the decaying baryon [P B = 0) we recover 
in the limit of a point-like meson * M(yf (£2 ~~ £2)) — > (27r) 3/2 ^/|(5(^2' - £ 2 ) the 
operational structure of the EEM (see eqs. |10] and [11]) , say 



(B'M\T\B) - 



^-L_5(P B - P B , - fc)(-73V37r)($ B /$ M |$ B * pair ) 



* B (li,6) (28) 



Hence it is obvious that we reproduce the EEM amplitudes by making in eq. 
( |TTD the replacements: 



2m n 



3i/, 



(2^)V2 m ^ 



7 



3^ 



(29) 
(30) 



In other words the compositeness of the meson reflects, with regard to the 
EEM, not only through the meson wave function but also in the transition ope- 
rator. 
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For the pion wave function we shall make use of a gaussian form fitted (Ra = 8 
GeV -2 ) to reproduce the root mean square radius of the pion [16| : 



^ {rl = ^^ eXP {-2W A ) (31) 

With an orientative purpose we compare it with the wave function extracted 
from the knowledge of the electromagnetic pion from factor for which (at low Q 2 ) 
we can very approximately assume vector meson (p-meson) dominance |J 



2 

in 

FAQ 2 ) = -rfai (32) 

The Fourier transform of F n (Q 2 ) gives the charge density from which the pion 
wave function can be derived: 

^(r) = Foo(r)^exp(-^) (33) 

We should keep in mind however that this expression is meaningful only at 
long distances 

In fig. 8 we draw the two normalized pion wave functions where the very 
different short-distance behaviour can be appreciated. 



4.1 The normalized NNtt form factor 

To get the normalized form factor defined previously in eq. (p0[), we compare 
the matrix element for p f— ■> p f 7r° calculated from i^Bar (eq. [14]) with the one 
provided by the QPCM, (eq. ^). The results appear in fig. 9. According to 
our expectations the inclusion of the pion structure seems to play the same role 
that the increasing of the baryonic size in the EEM model. In this sense half of 
the difference between the theoretical and the phenomenological results has been 
corrected for Vjj and Vm (for Vj the quantitative change is not very relevant). 

4.2 Fitting of 7 

Although the formal replacement in (|30|) establishes a relationship between f qq7C 
and 7 it should be handled with care since it is only valid in the point-like limit. 
If we extract 7 from /^^(O) in the same way as we did to get f gqn we have: 

i WO) f 1 + ^\ = « /(0) (34) 



(2m n y/ 2 m n V 2m N J 12tt 3 
where 

1(0) = J di di' di 2 di' 2 cos 2 sin 2 ^'2-esinV £ 5 £' 6 (l 



£ sin<; 
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* w (a) 



-1*^(0*^(0 - ^ | 3£cos0 2 cos0 + £'cos0 



\ 



'£sin<; 



COS# 2 0080*^,(^*^(0 - COS #2' 



£sin<; 



(35) 



is the Heaviside step function and the argument of the pion wave function is 



£ 2 cosS 



e 2 -2ecos0vr 2 -e 



sin 



(36) 



This gives for 7 the values shown in table 4. We observe for Vj a quite differ- 
ent value of 7 than for Vjj and Vni- This comes about due to the introduction of 
the pion structure; in the EEM model the constancy of f qq7T had to do with the 
fact that the baryonic quark structure contributes through model independent 
normalization integrals whereas in the QPCM case the presence of the pion wave 
function weights in a different form the quark structures giving provoking a model 
dependence in 7 that reflects the very different quark-model baryon wave func- 
tions, with and without three-quark force. Concerning the pion wave function 
its very short distance behaviour seems to affect the results very little. It turns 
out that this will also be the case for strong pion decays. For this reason we shall 
restrict hereforth our presentation to the gaussian wave function eq. ([H]). Let us 
notice that the values of 7 available in the literature are much lower than ours, 
the difference coming from the quark model used (harmonic oscillator) and from 



having fitted 7 in a different way. In ref. 



the vertex NNn was calculated 
analytically with the harmonic oscillator model. By choosing a value of R\ = 8 
GeV -2 for the pion wave function and a nucleon r.m.s. R 2 N = 6 GeV -2 (very 
close to the one obtained with Vj) the constant 7 would be 6.29 (6.97 with Vj). 



4.3 Strong Pion Decays 

The widths obtained with the 3 Pq model appear in table 5. We see that the 
introduction of the pion structure does not seem to represent any improvement 
of the results but rather on the contrary a worsening of the fit. 

As the specific form of the pion wave function is not very relevant we center 
our attention in the transition operator. Possible corrections to it may have a 
relativistic origin and can be enforced on the basis of general principles. First 
there is a normalization factor ./— which is related to the boost of the pion 
from its rest frame to the rest frame of the initial baryon. On the other hand, 
if one desires the quark momentum term to have the same properties as the k 
dependent one under exchanging the role of initial and final states it has to be 
associated to the difference of the energies of the initial and final states, which is 
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nothing but u w . In the non-relativistic limit, this factor tends to 2m q and it is 
therefore natural to introduce the factor ■^ JL - in front of the momentum term. 
This considerations can be implemented in eq. ( ]2"7j ) by introducing a factor 
and by multiplying the momentum (p^ 2 + p^) dependent term by 

Actually by making the additional replacement of the k coefficient | by (l + 
these changes are equivalent to modify the QPCM to get as the point pion limit 
the EEM. Adopting this point of view we define our 3 Po modified QPCM through 
the operator: 



o mod (g,6) = sj—yr 




2m q V 3 





*J,W§(6 - 6)) exp I <y - 6) 

ex P I m/hT^S ~ 6) ) exp [ i\ -k& I (37) 



By repeating the calculational process as explained in section 4.1 and 4.2 
the improvement of the fit with this modified QPCM is spectacular as seen in 
table 6 (the normalized form factor hardly changes and we do not draw it again). 
Except for A(1600) — > Ntc and A^(1535) — > Aw where the strong cancellation 
between derivative and non-derivative terms makes the result very sensitive to 
their precise values all the widths are reproduced withing a factor 2. (Let us 
remark that the experimental interval has been taken as the most-restrictive 
one, i.e. by applying the corresponding variable decay percentage to the central 
value of the width). 

The results for the transitions involving the nucleon (or A particle) and its 
radial excitations, A^(1440) (or A(1600)) should be noticed. The finite size of 
the pion makes them sizeable and in rough agreement with experiment. This is 
achieved without relying on (^-) 2 corrections which in the EEM would probably 
allow one to get a similar result. Indeed, the orthogonality of wave functions, 
which explains the low transition rates given in table 3, does not apply anymore 
when the operator nature of these corrections is accounted for. 

The overall results for the three potentials do not differ very much. Added 
to the NNk form factor prediction and keeping in mind that other corrections, 
as the ones coming from the A-width in the final states, could be evaluated, the 
small remaining discrepancies with experimental data might then have to do with 
a fine tuning improvement of the quark potential on the base of increasing the 
quark model radii. 
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5 Summary 



A refinement of a non-relativistic quark model that by means of the incorporation 
of a three quark potential is able to reproduce the nucleon and delta low energy 
spectra has been carried out avoiding some conceptual difficulties associated to 
the very singular form of the three-body force and to the values of the potential 
parameters in the original version. The predictions for static properties point 
out the need of incorporating (qqq qq) components into the wave function. An 
effective manner to do this when studying strong decay processes is through the 
transition operator defining the calculation scheme either via a modification of 
the coupling constant as in a elementary emission model (EEM) or additionally 
via the explicit implementation of qq pair creation as in the quark pair creation 
models (QPCM). 

We do find some sensitivity to the way the 7r-decay process occurs at the 
microscopic level. So the 3 P QPCM largely overestimates majority of the decay 
rates. Discrepancies were found to have an origin in relativistic effects. Starting 
from this observation, the QPCM was modified by incorporating those corrections 
which stem from general principles, making it closer to the EEM. The quality 
of the results we thus obtain seems to support our prescription. Despite large 
differences in the short distance description of baryons used in our calculations, we 
did not find any reason to discriminate them by looking at the n decay properties 
once the experimental kinematics is considered. Introducing a three body force in 
the description of baryons makes the agreement with experiment slightly better in 
some cases, slightly worse in others. Probably, the observable we examined is not 
appropriate to extract any correlation between the predictions for the spectrum 
and the decay widths. Anyhow as a general conclusion it can be established 
that a unified description of the spectrum and the strong baryonic decays in a 
'non-relativistic' scheme seems plausible. 

This work has been partially supported by CICYT under grant AEN93-0234 
and by DGICYT under grant PB91-0119-C02-01. F.C. acknowledges the Min- 
isterio de Education y Ciencia for a FPI fellowship. B. D. has been partially 
supported by Conselleria d'Educacio y Ciencia of the Generalitat Valenciana. 

Appendix A 

All we can know about the baryonic structure is contained in the intrinsic wave 
function ^(^1,^2) depending on the six independents coordinates £1,^2- Alterna- 
tively one can define hyperspherical coordinates in the dimension 6-hyperspace 
as an hyperradius £ defined through 
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£ 2 = £i 2 + £ 2 2 (A.l) 

an a set Q of 5 angles: £1 (the spherical angles of £1), £2 (the spherical angles of 
£2), and given by 



6 

6 



£ sin ( 

£ COS( 



(A.2) 



such that <i£ 2 = <i£ £ 5 d£l = d£ £ 5 c?0 sin 2 </> cos 2 0<i£i d^2 

In terms of Q one can construct a complete set of basis functions on the 
hypersphere unity. There are the hyperspherical harmonics (HH) Yrjn(fi), char- 
acterized by 5 quantum numbers denoted by [K] = [K, li, mi, I2, m 2 ] ||. 

Then the spatial part of the intrinsic wave function can be expanded as: 



satisfying the Schrodinger equation: 



[K] 



(A.3) 



(A.4) 



where H int stands for the sum of the internal kinetic and the potential energy. 
Explicitly by defining reduced radial wave functions as 



and defining the wave function vector 

one has the set of coupled Schrodinger equations: 



(A.5) 



(A.6) 



1 



d 2 



v K 



where V is the matrix of the potential 



and vk is the diagonal matrix 



■1.1 



(A.7) 



(A. 



(A.9) 



The set of eqs. (|A.7 ) can be solved numerically by Numerov integration, 
turning out to be convenient the introduction of the dimensionless variable: 
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he 

in terms of which we can finally write 



x = 4^ (A.10) 



d 2 v K V(x) E] . , , 4 . 

+ — -- XW=0 (A.ll) 



dx 2 x 2 



m m 



Let us note that under a variation of the quark mass m — > m! (x — > x' = — x) 
( |A.11| ) transforms into: 



d 2 V(x') _ E/_ 

x' 2 ml m' 



where 



— 771 

V(x') = —V(x) (A.13) 
ml 

Tfl 

E' = —E (A.14) 

X(x') = X (x) (A.15) 



So, by redefining the parameters of the potential to satisfy (|A.13 ) we can 



obtain an energy spectrum rescaled by a factor without any change in the 
wave functions. 

For physical purposes it is more convenient to use linear combinations of HH 
with total orbital angular momentum L and definite symmetry, ^[x^Lnetiy] (^) ^ n 
terms of which the total intrinsic wave function reads: 



*b(£i,6)$bOS,M5;I,M7 



JJz 



J2 {L l S l J\M l M St J z 

[symmetry] 



^;SL ry] $[— y](Si, M s - I t , M h )] (A.16) 



having to be totally symmetric as required by the symmetrization postulate. 
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Table captions 

Table 1. Fitted values of the parameters of the two-body and three-body 
potentials. 

Table 2. Magnetic moments (in nuclear magnetons) and radii (if fm 2 ) of the 
nucleon. 

Table 3. Decay widths (in MeV) with an EEM model. Experimental data from 

0- 

Table 4. Values of 7 obtained from the NNw vertex. 

Table 5. Decay widths (in MeV) with a 3 Po model. Pion wave function as 
given by eq. fl3T| ) and 7 taken from table 4. 

Table 6. Decay widths (in MeV) with a modified 3 P QPCM. rion wave 
function given by eq. (|3l"l). 
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Figure captions 



Figure 1. Relative energy spectrum for positive parity nucleon and delta and 
negative parity nucleon states. Solid lines correspond to the predictions of 
Vj and the shaded region whose size represents the experimental width to 
the experimental data J7|. 

Figure 2. Idem fig. 1 with Vi — > Vn. 

Figure 3. Idem fig. 1 with Vj — > Via- 

Figure 4. Spatial symmetric component of the nucleon reduced wave function 
The solid line corresponds to Vj, the dashed-dot line to Vn and the 
dashed line to Vni- 

Figure 5. A baryon B decays through pion emission by one of its quarks. 

Figure 6. Normalized NNir form factor. The dot line corresponds to the phe- 
nomenological parameterization of ref . [15j . The solid line, dashed-dot and 



dashed line correspond to the predictions of V/, Vn and Vni respectively. 

Figure 7. A quark and antiquark from the vacuum recombine with the quarks 
of the initial baryon to give rise to the final baryon and pion. 

Figure 8. Normalized pion wave function. The solid line is obtained from the 
Fourier transform of the pion form factor (^) The dashed line corresponds 
to the gaussian form eq. (PTT). 



Figure 9. Normalized NNir form factor in the 3 Pq model. The line corresponds 
to the phenomenological parameterization of ref. [15|. For the pion wave 
function, eq. (gj) with R\ = 8 GeV" 2 . 
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Vi 


Vii 


t r 

VIII 


m u = m d (GeV) 


0.337 


0.355 


0.320 


y(COUL) 


k (GeV fm) 


0.1027 


0.289 


0.321 




K a (GeV fm) 
r (fm) 


0.1027 
0.4545 


0.049 
0.40 


0.044 
0.49 


y(CONF) 


a 2 (GeV- 1 fm) 


1.063 


4.570 


4.124 




V (GeV" 2 fm" 6 ) 
m (GeV) 




-61.63 
0.25 




T/(3) 

V II 


Vo (GeV) 
A (fm) 






-35.5 
0.25 



Table 1 
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Vi Vu V in 



Exp. 



fi p 2.76 2.64 2.93 2.79 

fi n -1.83 -1.76 -1.95 -1.91 

(r 2 N } m 0.218 0.128 0.115 

(r2) Ch 0.238 0.133 0.112 0.74 ± 0.02 

(Oc -0.02 -0.005 -0.004 -0.119 ± 0.004 



Table 2 



23 



Vi V n V in Exp. 

A(1232) 
7V(1440) 
7V(1440) 
A(1600) 
A(1600) 
7V(1520) 
7V(1520) 
7V(1535) 
iV(1535) 



Table 3 



Vi 


Vn 


Vin 


7 7.02 


9.78 


11.04 



Table 4 



-> Ntt 79.6 72.1 

Nit 3.4 0.17 

^Att 7.1 17.6 

-> Nn 20.1 94.1 

^Att 2.85 0.10 

-> A% 61.8 22.3 

-> Att 78.0 56.1 

-> A% 240 149 

-> Att 9.7 8.3 



67.7 115-125 
0.01 210-245 
24.1 70-105 
148.6 35-88 
0.08 140-245 

17.8 60-72 
55.3 18-30 
117 53-83 
9.1 <1.5 



Vi V n Vjjj Exp. 



A(1232) -> Ntt 167 210 241 115-125 

iV(1440) -> iVvr 452 1076 1509 210-245 

7V(1440) -> Att 66.5 228 337 70-105 

A(1600) -> Nn 19.8 0.53 7.0 35-88 

A(1600) -> Att 255 498 656 140-245 

iV(1520) -> iVTr 268 319 353 60-72 

iV(1520) -> Att 532 999 1203 18-30 

7V(1535) -> JVtt 429 464 410 53-83 

iV(1535)^A7r 28.1 74.0 107 <1.5 



Table 5 
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Vj V n V IH Exp. 



A(1232) -> Nit 88.6 

iV(1440) ^ iVvr 114 

7V(1440) -> Att 27.6 

A(1600) -> iVTr 2.1 

A(1600) -> Att 62.0 

7V(1520) -> iVvr 95.1 

iV(1520) -> Att 45.9 

7V(1535) -> iVvr 49.2 

iV(1535) -> Att 15.3 



112 123 115-125 

307 469 210-245 

116 183 70-105 

2.5 5.6 35-88 

121 177 140-245 

105 118 60-72 

75.2 111 18-30 
44.1 55.8 53-83 

36.3 52.1 <1.5 



Table 6 
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